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We present measurements of tunneling magneto-resistance (TMR) in single-wall carbon nanotubes 
attached to ferromagnetic contacts in the Coulomb blockade regime. Strong variations of the TMR 
with gate voltage over a range of four conductance resonances, including a peculiar double-dip 
signature, are observed. The data is compared to calculations in the ’’dressed second order” (DSO) 
framework. In this non-perturbative theory, conductance peak positions and linewidths are affected 
by charge fluctuations incorporating the properties of the carbon nanotube quantum dot and the 
ferromagnetic leads. The theory is able to qualitatively reproduce the experimental data. 
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I. INTRODUCTION 

Controlling electronic spin in nano-scale circuits is 
a long-lasting challenge on the way to fast-switching, 
energy-efficient building blocks for electronic devices. 
To this end, spin-dependent transport properties have 
been investigated in a wealth of low dimensional sys¬ 
tems, e.g., mesoscopic magnetic islands if], 2DEGs [2], 
InAs nanowires [3], graphene 0] and fullerenes [5]. Car¬ 
bon nanotubes (CNTs), being thin, durable and high- 
throughput wiring, allow coherent transport of electronic 
charge and spin and are promising candidates for future 
spintronics applications fEy. While control and scalabil¬ 
ity of CNT-based nanocircuits still pose significant chal¬ 
lenges, devices where single carbon nanotubes (CNTs) 
are contacted to ferromagnetic leads can be produced 
with standard lithography methods: spin valve exper¬ 
iments were performed on single-wall [7HIP] (SWCNT) 
and multi-wall [11H15] carbon nanotubes in various elec¬ 
tron transport regimes. In most cases, a spatially con¬ 
fined quantum dot is coupled to ferromagnetic electrodes. 
Electronic transport across CNT quantum dots can take 
place in different regimes: Depending on the relative 
magnitude of coupling strength, temperature and charg¬ 
ing energy, this ranges from an opaque Coulomb-blockade 
regime |16H19I . to an intermediate coupling regime with 
lead induced energy level shifts [201422] , to a strongly cor¬ 
related Kondo regime [25H2I)] . For highly transparent 
contacts, in contrast, the dot behaves essentially like an 

electronic wave guide EH 123- 

in our work, we focus on the conductance of a car¬ 
bon nanotube quantum dot weakly coupled to ferromag¬ 
netic contact electrodes, recorded for parallel (G p ) and 
anti-parallel (G ap ) contact magnetization, respectively. 
G p and G ap define the so-called tunneling magneto- 
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resistance (TMR) gZHMj: TMR = (G p /G ap ) - 1. Ex¬ 
perimentally, the TMR has been shown to be strongly 
gate dependent mm- We report on shifting and broad¬ 
ening of conductance peaks resulting in specific dip-peak 
and dip-dip sequences in the TMR gate dependence. Our 
data covers a range of four Coulomb resonances with ex¬ 
tremal TMR values of —20% to +180%. 

The pronounced resonant structure of the conduc¬ 
tances G p and G ap leads to large TMR values if the 
positions and widths of the resonances depend on the 
magnetization configurations p and ap. Thus, various 
mechanisms have been proposed which induce a shift of 
the energy levels of the quantum dot, and thus the of 
the resonance peaks, depending on the magnetization of 
the contacts. Those are spin-dependent interfacial phase 
shifts [7] or virtual charge fluctuation processes mm- 
The effect of spin polarized leads on the resonance width 
have been described in j31| for a resonant single level 
junction. Interestingly, a negative TMR is predicted for 
asymmetric couplings to the leads. An attempt to ac¬ 
count for broadening in the presence of Coulomb inter¬ 
actions was discussed within a self-consistent approach 
based on the equation of motion (EOM) technique |2~0| . 
The EOM was applied to model TMR data reported for 
a SWCNT [7] for a model with spin-dependent interfacial 
phase shifts. 

Here we discuss a transport theory which naturally in¬ 
corporates the effects of spin polarized leads on the po¬ 
sition and width of conductance resonances in the pres¬ 
ence of strong Coulomb interactions. It is an extension 
of the so-called dressed second order (DSO) transport 
theory, recently developed for normal leads [22], to the 
case of spin-polarized contacts. This theory accounts for 
energy renormalization and broadening of the peaks in 
linear conductance due to charge fluctuation processes. 
We show that the charge fluctuations also affect trans¬ 
port through excited states in the non-linear conductance 
regime. This observation is in agreement with previ¬ 
ous reports on tilted co-tunneling lines in CNT quantum 
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FIG. 1. (Color online) SEM picture of a chip structure simi¬ 
lar to that of the measured device. A carbon nanotube on a 
positively doped silicon substrate capped with 500 nm SiC >2 is 
contacted by two Permalloy stripes, one of which is exchange- 
biased by a FeMn layer. On top, the stripes are protected by 
palladium. Gold is used for the bond pads and the connec¬ 
tions to the nanotube contacts. 


FIG. 2. (Color online) Differential conductance versus bias 
and gate voltage of a selected region measured at 300 mK and 
B = 0. The numbers in the Coulomb blockade regions denote 
the number of electrons in shell n on the quantum dot. Arrows 
indicate the first excited state crossing the source (left) and 
drain (right) lines in the vicinity of the state with one extra 
electron (N = 4n + 1). 


dots [32]. A qualitative agreement with the experimental 
findings is obtained. 

This paper is structured as follows. We first present 
the measurement details and experimental data in Sec. [TT] 
In Sec. m we introduce the so-called dressed second or¬ 
der theory (DSO) (22! in the reduced density matrix 
transport framework and address its implications on non¬ 
linear conductance and TMR. Finally, in Sec. |IV| we pro¬ 
vide a comparison between experimental data and results 
from the DSO and draw our conclusions in Sec. El 

II. EXPERIMENT 
A. Sample preparation 

For the purpose of measuring TMR in CNTs, one needs 
to interface the nanotube to two ferromagnetic contacts 
with a different switching field. The conductance, being 
sensitive to the magnetization in the leads, changes when 
the polarization of one of the contacts is reversed by an 
external magnetic field. It has been shown that NiFe 
is well suited as a material for the electrodes of CNT 
spin-valves [33] : the alloy shows a distinct switching be¬ 
havior as a function of the applied magnetic field and 
the interface transparency between NiFe and the CNT is 
comparable to that of Pd. The structure of one of the 
devices we realized for this purpose is shown in Fig. |Tj 
On an oxidized silicon substrate (500nm Si02) a carbon 
nanotube is grown by chemical vapor deposition. The 
nanotube is located by atomic force microscopy and two 
NiFe (80:20) leads, 20 nm in thickness, are deposited at 
a distance of 1 /im on top of the nanotube by sputtering. 
On one of the two contacts, 40 nm of anti-ferromagnetic 
FeMn (50:50) is sputtered to bias the magnetization of 


the underlying NiFe contact. The hysteresis loop of this 
contact is expected to be shifted with respect to the pure 
NiFe contact by virtue of the exchange bias effect !331 . 
A 20 nm protective layer (Pd) covers the leads from the 
top. The switching of the exchange biased contacts was 
confirmed independently prior to the measurement using 
SQUID and vibrating sample magnetometer techniques. 


B. Measurement 

An electronic characterization of the quantum dot 
at 300 mK and at zero magnetic held shows a regular 
Coulomb blockade behavior (Fig. [2]). The data yield a 
gate conversion factor a = 0.29 and a charging energy 
of E c = 6.1 meV (see Eq. 0 ). The sample does not 
exhibit a clear four-fold symmetry in the peak height or 
peak spacing as expected for a carbon nanotube quantum 
dot. Consequently, we are not able to label the Coulomb 
blockade regions with a value of the electronic shell fill¬ 
ing n in a definite way. The assignment of the number 
of electrons to the experimental data in Fig. [2] is done in 
agreement with the theoretical predictions in Sec. |III E 

Having a closer look at Fig. [2] we can identify an ex¬ 
cited state transition at 1.4 meV parallel to the source 
line (left arrow) and at ~ 1.8 meV parallel to the drain 
line (right arrow). The energy scale of this excitation 
stays approximately constant over a range of at least 
six resonances, as can be seen from measurements over 
a broader gate range. The quantization energy e(n) of 
a CNT shell n is a direct consequence of the electron 
confinement along the nanotube. It yields a mean level 
spacing eo = e(n + 1) — e(n) oc Hvf/ttL, where L is the 
CNT length. It is thus reasonable to identify the first 
excitation with the confinement energy eo equivalent to 
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a lateral confinement of 1.1 /xm for a Fermi velocity of 
800km/s [35], a value close to the contact spacing of 
1 /jin. The asymmetry of the line spacing with respect to 
source and drain suggests a gate-dependent renormaliza¬ 
tion [32] of the CNT many-body addition energies in the 
presence of ferromagnetic contacts. We show in Sec. m 
that this can be a direct consequence of charge fluctua¬ 
tions in the presence of contact magnetization. 

Electron transport measurements at 300 mK show a 
significant switching behavior. In Fig. [3] the conductance 
across the CNT quantum dot is plotted against the mag¬ 
netic field directed parallel to the stripes, i.e., along their 
easy axis, as indicated in the inset to the figure. The steps 
in the signal can be interpreted as the magnetization re¬ 
versal of the contacts, as sketched in the figure. Sweeping 
the magnetic field from negative (—100 mT) to positive 
values, one of the contacts switches at H = H s u , result¬ 
ing in a configuration with anti-parallel polarization of 
the majority spins of the two contacts. This results in a 
drop of the conductance signal. Upon increasing the field 
further, the second contact is supposed to switch and the 
conductance should recover. The second switching event 
was not observed in the present sample. Sweeping back 
from positive to negative field, the conductance recovers 
at Hg^. The two values H s d / u characterize a hysteresis 
loop with a coercive field H c = H s u — H s <3 and an ex¬ 
change bias H ex = (U s ,u + H Sy d)/2. At B = 0 the two 
contacts are always in a parallel configuration, because 
the coercive field of the switching contact is smaller than 
the exchange bias. 

Measurements of the conductance performed at zero 
magnetic field require Ar ast ~ 100 ms per data point 
and will be called the fast measurements in the follow¬ 
ing. Contrarily, in slow measurements, each conductance 
data point is obtained from magnetic field sweeps with a 
duration of Af slow ~ 20 minutes at constant gate voltage 
(compare Fig. [ 3 ]). We then identify H s from a step in the 
conductance signal and take the average over 100 points 
on either side of the step to extract the conductance in 
the parallel and anti-parallel configuration, respectively. 
This is repeated for 250 values of the backgate potential 
in the range between 8.126 V and 8.201V. In Fig. [4] the 
TMR as a function of gate voltage is shown together with 
the conductance at parallel contact polarization. In this 
slow measurement, we obtain conductance peaks with 
a height of 0.15e 2 //i and a full width at half maximum 
(FWHM) of r ~ 0.7 meV. Comparing these values to a 
height of 0.3 e 2 /h and a width of 0.4 meV obtained from 
the fast measurement at B = 0 we conclude that the peak 
conductance in the data from the slow measurement is 
substantially suppressed. We will discuss this deviation 
in Sec. [IV] It is remarkable that besides huge positive 
(180%) TMR values, negative regions occur prior to the 
peak in the TMR curve in the first two resonances while 
for the last two the value drops again, forming two dips 
in sequence. Again this will be discussed in more detail 
in Sec. ITVl 



FIG. 3. (Color online) Differential conductance plotted ver¬ 
sus magnetic field at V g = 8.1737 V, Vb = 0 and 300 mK. The 
solid red curve was recorded with increasing field, the dashed 
blue curve with decreasing field. Small pictograms indicate 
possible orientations of the majority spins in the contacts. 
The switching of one of the two contacts at -ff s , u /d is high¬ 
lighted with arrows at the bottom for both sweep directions. 
The coercive field is indicated by H c and the exchange bias by 
H ex . Inset: Orientation of the external field B with respect 
to the CNT and the leads. 
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gate voltage (V) 


FIG. 4. (Color online) Differential conductance and TMR as 
a function of gate voltage measured over four resonances (slow 
measurement, see text). The conductance is measured at par¬ 
allel polarization of the contacts. The TMR graph shows a 
dip-peak sequence over the first two resonances and a quali¬ 
tatively different double-dip feature at the last two. 


III. THEORETICAL MODELING 

We proceed by presenting a theoretical framework ca¬ 
pable to reproduce the transport data from the previous 
section. In particular, the connection between the theory 
and the resulting shape of the TMR curve will be dis¬ 
cussed in detail. In order to be able to account for a gate 
dependence of the TMR, the transport theory should be 
able to incorporate the influence of the ferromagnetically 
polarized leads on the positions of the linear conductance 
maxima as well as on the width of the conductance peaks. 

Noticeably, the commonly used perturbative descrip¬ 
tion of the Coulomb resonances predicts temperature 
broadened peaks and maxima whose positions are solely 
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is the Hamiltonian of an ensemble of non-interacting elec¬ 
trons in the leads l = s/d with wave vector k and spin er. 
The operator Ci a k (cJ Tak ) annihilates (creates) an elec¬ 
tron with energy e; CT k- The second part, 


- 1 ~ 2 

H d = -E C N 


E 


e(n ) + tct 


A S o (n) 


N„ 


- eaV g N + Hgxt A » 


(1) 


FIG. 5. (Color online) Left panels: Schematic drawing of the 
lead induced, polarization dependent, modification of posi¬ 
tion (a) and width (b) of a peak in the conductance across a 
quantum dot as a function of the gate voltage. Right panels: 
As a consequence of the level shift (a) and level broadening 
(b), the corresponding TMR signal exhibits a characteristic 
dip-peak (c), or dip-peak-dip (d) feature. 


determined by the isolated quantum dot spectrum im¬ 
plying a constant, positive TMR [2l] . 

A transport theory accounting for charge fluctuations 
non-perturbatively was shown to shift the quantum dot 
energy levels depending on the magnetization configu¬ 
ration of the leads m- The qualitative effect of the 
renormalization is depicted in Fig. |5](a) : The peak in the 
conductance G p in presence of leads with parallel spin 
polarization is shifted with respect to the one in G ap , the 
conductance in the anti-parallel case. This shift yields a 
characteristic dip-peak feature in the TMR signal, similar 
to what was observed in Ref. [7[ Yet, this theory cannot 
account for the double-dip like TMR signatures visible in 
our data (see Fig.EJ V g ~ 8.19 V and V g ~ 8.17 V). These 
require additionally a change of the resonance line-width 
when switching from the parallel to the anti-parallel con¬ 
figuration, as shown in Fig. [5](b) and also observed lately 

in Ref. [501 . 

In the following, we discuss how to theoretically ac¬ 
count for broadening and renormalization effects, to low¬ 
est order in the coupling T, within the recently proposed 
“dressed second order approximation” (DSO). The DSO 
has been discussed in Ref. [22] for the single impurity An¬ 
derson model with normal metal leads, where it has been 
shown to correctly capture the crossover from thermally 
broadened to tunneling broadened conductance peaks. 
Here we present its generalization to a multilevel system 
coupled to ferromagnetic leads. 


A. Hamiltonian 

We treat the system as an isolated quantum dot cou¬ 
pled to metallic leads. The Hamiltonian of such a system 
reads H = Hr + Hd + Ht- Here, 

Hr = Gr7k:cj tT kGo‘k 
Icr k 


describes the electrons on the CNT quantum dot in terms 
of the quantum numbers n (shell), spin a and valley r. 
Here we used N nrcr = d| lro .(i nrcr , with the fermionic dot 
operator d nTa and N = ^ nT0 . N„ rcr , the total dot oc¬ 
cupation. For our purposes, it is sufficient to account 
for Coulomb interaction effects in terms of a capaci¬ 
tive charging energy E c . Short range exchange contri¬ 
butions are neglected here. The symbols r and a repre¬ 
sent the eigenvalues ±1 of the states with quantum num¬ 
bers K,K' andt,i, respectively. In the CNT, a non-zero 
spin-orbit coupling Aso can lead to the formation of de¬ 
generate Kramer pairs |26| . Notice that, for simplicity, a 
valley mixing contribution is not included in Eq. ([I]) , as it 
would not affect the main conclusions drawn in this work. 
Hence, the valley degree of freedom is a good quantum 
number to classify the CNT’s states [30]. The next to 
last part of the Hamiltonian Hd models the effect of an 
electrostatic gate voltage V g scaled by a conversion fac- 

^ p/a 

tor a. Finally, H oxt accounts for external influences on 
the dot potential, e.g., stray fields from the contacts and 
the external magnetic field used to switch the contact 
polarization. 


The ground states of shell n have An + a (0 < a < 3) 
electrons and will in the following be characterized by the 
quantum numbers of the excess electrons with respect to 
the highest filled shell n — 1. For instance, the quan¬ 
tum dot state labeled by \K^\n) contains 4 n electrons 
plus one additional electron in the (K, f) state. Includ¬ 
ing states with 4n — 1 and 4n + 5 electrons we end up with 
6 ground states with different degeneracies (see Tab. IJ 
left column). In total we consider a Fock space of di¬ 
mension 24 if the four-fold degeneracy is not lifted by 
a sufficiently large spin-orbit coupling Aso- The extra 
states with occupation An — 1 and An + 5 are included 
to allow for charge fluctuations in and out of the shell n 
under consideration. Conversely, for large enough spin- 
orbit coupling the dimension of the Fock space is reduced 
to 10, see Tab. [TJ right column. Judging from the stabil¬ 
ity diagram in Fig. [2] and from data over a greater gate 
range where we see no two-fold pattern in the spacing 
of the excited state lines, we consider the configuration 
on the left side in Tab. [T| to be more likely. For a com¬ 
pact notation, the shell number will in the following be 
neglected from the state ket if not necessary. 


Quantum dot and metallic leads are coupled perturba- 
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We 1 Aso < max{k B T, 70 } 

Aso > max{k B T, 70 } 

-1 


|A't,A4,A"t;n-l) 


\Kt,Kl,K'l;n-l) 

1 ) 

\Kl,K'^,K'l,n-l) 

\Kt,Kl,K'l,n-l) 

0 

1 n) 

\n) 

1 

|A't;n) | Ki-n) 

\K' t;n> A"4_; n) 

\Kt,n) \Kl-n) 

2 

| ATT, A't; n) 
\Kt,K'i-n) |Af,A"t;n) 

1 Ki,K'i-n) \K'tK'i-n) 

m ,Kl,n) 

3 

\Kt,Ki,K't,n) 

| At, Ki, A't; n) 


|A't,Af,A't;n) 

|A't, A't, AT'4,; n) 
\Ki,K't,K'i-n) 

|At,A4,A'f;n) 

4 

\n + 1> 

|n+ 1) 

5 

|A't; n + 1) [A'4_;n+ 1) 
|A"t;n+l) \K'l-n+l) 

\Kfrn + l) [A'Rn + l) 


TABLE I. The set of allowed electronic ground states C of 
the CNT with N electrons for large (right) and small (left) 
spin-orbit coupling Aso- The degeneracy of the configuration 
depends on the magnitude of Aso- In the first column, the 
excess electron number Wei = N — An is reported with respect 
to the number 4n of electrons in the filled (n — l)-th shell. 


tively by a tunneling Hamiltonian 


given in second order by 


K\ 


bb 


= Vrf, =y 1 lim 

l,ba z_^ A _ >0+ 


, TfWW 

E b a -e + iX 


h. c., 


where F^ ba is the corresponding tunneling rate. The 
function /f(e) with p = ± is defined as /,(e) = 
[1 + exp{±/?(e — pi)}]- 1 , where /3 is the inverse tem¬ 
perature and pi the lead’s chemical potential. Hence, 
/, + (e) = fi(e) is the Fermi function and describes the 
occupation probability in lead l. In general, p = ±1 if 
the final state | b) has one electron more/less than the 
initial state |a). The energy difference between final and 
initial dot configuration is given by = E b — E a = 
E b - E a - eaV g (N b - N a ). Finally, 

7?°( e ) = 7io-(6,a)( e ) = \Ti\ 2 'Di cr (e) 

is a spin-dependent linewidth defined in terms of the tun¬ 
neling amplitude T) and of the spin-dependent density of 
states 2?zcr(e). A Lorentzian provides a cut-off for the den¬ 
sity of states at a bandwidth W. The notation a(a, b) 
indicates that the spin a of the electron tunneling out 
of/onto lead l depends on the spin configuration of the 
initial state a and the final state b of the quantum dot. 
It is convenient to introduce the spin-resolved density of 
states of lead l at the Fermi energy 


Ht — ^ ' Tikno-rdncn-Cikcr T h. C., (2) 

ZkntTT 

with a tunnel coupling T/kno-T generally dependent on the 
quantum numbers of both leads and quantum dot. In the 
following, for simplicity, we assume that Ti^ naT = T). 


B. The reduced density matrix within the dressed 
second order (DSO) approximation 

We describe the state of our system by the reduced 
density matrix p = Tra{/5 tot }, obtained by tracing over 
the possible configurations of states in the reservoirs, as¬ 
suming that they are in thermal equilibrium. For the 
quantum dot itself we suppose that it reaches a steady 
state characterized by p = 0. The corresponding station¬ 
ary Liouville equation reads [22j 

0 = - iy SabSa’V (Ea - E' a )p aa , + ^ K$p aa ', (3) 

aa' aa r 

in terms of matrix elements p ab = (a|/S|6) of p in the 
eigenbasis of the quantum dot. The superoperator K 
connects initial states | a), |a') to final states |6) and | b') 
at a certain order in the perturbation Ht- 

The calculation of the kernel elements is performed 
along the lines of Ref. [22] As an example, the element 
connecting the states |6), \b') = |&) and |a), |a') = |a) is 


V la =V l(T (e F )=D 0 (l + aP l )/2 (4) 

where Pi = (Wf — £> 4 )/( T >4 + V 4 ) is the polarization of 
lead l. The couplings |T )| 2 we define in the same spirit as 

iy d | 2 = |Io| 2 (l±a)/2, (5) 

using the parameter a to tune the asymmetry in the cou¬ 
pling to the leads. We will in the following use the fac¬ 
torization 

7i<r(b,a)(£F) = 70*4(7, (6) 

where we collect the lead and spin independent prefactors 
in an overall coupling strength y 0 = -D 0 |T 0 | 2 and include 
the dependence on spin and lead index in the dimension¬ 
less parameter ni a , where K i° = 1- Note that 70 is 
related to the level broadening r 0 by T 0 = 27r7o. 

In Fig.^a), a diagrammatic representation of one con¬ 
tribution to the second order kernel is shown for the case 
of | a) = |0) and | 6 ) = |rtr). The fermionic line connect¬ 
ing the lower to the upper contour carries indices e, a 
which fully characterize the nature of the electron tun¬ 
neling between lead l and quantum dot. The direction of 
the arrow further specifies if the electron tunnels out of 
(towards lower contour) or onto (towards upper contour) 
the dot. 

Beside this lowest (second) order contribution, we con¬ 
sider all diagrams of the structure shown in Fig. [ 6 ^b) . 
The selected diagrams contain arbitrary numbers of 
uncorrelated charge fluctuation processes (bubbles in 
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FIG. 6. (Color online) Diagrammatic representations of the 
contributions to the rate r+ TCT0 in second order (a), and an 
example of diagrams included in the DSO (b). In the latter 
case, the fermion line (blue) from the second order theory is 
“dressed” by charge fluctuation processes. The labels below 
the fermion lines denote energy and spin of the particle tun¬ 
neling from/onto the lead. Note that the diagram is read from 
right to left, i.e., the initial state |0) can be found on the right 
and the final state |rcr) on the left. 


for a state b that can be reached by an in-tunneling pro¬ 
cess from state a, and 

rC o6 =^:/de^( e )/r(e) (8) 


for an out-tunneling process b — > a. Note that we intro¬ 
duced a tunneling-like density of states (TDOS) 


,ba . 


:*) = 


7 j (e) Im(£“ a (e)) 


[Im(£ 6a (e))] 2 + [e - E b a + Re(£ ba (e))] 2 ' 


(9) 


We refer to the contribution £ ha in the denominator of 
the TDOS as a self energy that infers from the contribu¬ 
tions of all possible charge fluctuations connected to the 
initial, a, and final, b , states in the state space given in 
Tab. [T| Explicitly, 

£ ba (e)= E ( 10 ) 

cE{6,a} 

c'ecf 

with the sets C^, given by 


Fig.§. During the charge fluctuation, the dot state on 
the upper contour has one charge less or more compared 
to that of the final state \ra). Hence, the virtual state is 
either the state |0) or one of the many (see Tab. [T]) doubly 
occupied states. On the lower contour, the fluctuations 
take place with respect to the initial state |0). Exam¬ 
ples of charge fluctuations in the case of initial state |0) 
and final state |A'f) are shown in Fig. [?J Summing all 
diagrams of this type yields the DSO rates : 

r A = ^/< 

(7) 


C± a := {d : N c , = N b/a ± 1 A 4n - 1 < N & < 4n + 5}. 

( 11 ) 

The sets are shown in Fig. [?] for the states |a) = |0) and 

\b) = \Kt). 

The summand 


c' (b/a) 


%a 


(e) = E / dw 


ui) 

±pu + e- E%/* + iq 


accounts for a transition from b or a to a state d, with d £ 
C^,. Performing the integral, we arrive at an analytic 
expression for the contributions to the self energy, i.e., 


c '( b / a )/A _ ' 


%a 


(c) = E 7E /0) (e) {iTr/f (±P(E C ^, - e)) ± [^^(W) - Re [*(°> (i( W ±- e)))]] } , (12) 


where 4d 0 )(a;) = T ,< ^ o ^(0.5 + i/27rkBT) and is the 
digamma function. Note that the dependency on the 
bandwidth drops out due to the alternating sign of the 
contributions from the upper and lower contour in the 
summation in Eq. (101. Having calculated the self en¬ 
ergy, we are now able to collect all rates according to the 
transitions in our state space, and solve the stationary 
Eq. @ to obtain the occupation probabilities p aa = P a . 
Within the steady state limit we can neglect off-diagonal 
entries p ba if they are among non-degenerate states [21] . 
According to Tab. [I] the CNT spectrum can be spin and 


valley degenerate. However, the tunneling Hamiltonian 
([2]) conserves the spin during tunneling, and thus spin 
coherences are not present in the dynamics. Here, for 
simplicity, orbital coherences are neglected as well EH- 


C. Current within the DSO 


The current through the terminal l can be written in 
terms of the difference of in- and out-tunneling contribu- 
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tions at the junction |38j : 


/z(V b ) = 


e 

27 tH 


£ / de [Pa(V b )/+(e)-P c (V b )/r(e)] Vb). 

aec J 
ce C+ 

(13) 


where V b is the bias voltage applied between the two 
contacts, and C is the set of all possible configurations 
(see Tab. [I]). In general, the populations can be expressed 
in terms of rates via the Liouville equation ([3]) and a 
closed form for the current and, consequently, for the 
conductance can be found. This is straightforward if two 
states are connected by pairwise gain-loss relations [58] . 
For the case of the single impurity Anderson model, for 
example, a compact notation of the conductance can be 
given [22]. In this work, the conductance data from the 
model is calculated numerically [35]. 

The width of a resonance in conductance with respect 
to the gate potential is determined by the populations, 
the TDOS which has a form similar to a Lorentzian, and 
by the derivative of the Fermi functions. Note that the 
populations are themselves a function of the rates and 
therefore are also governed by the resonance conditions 
of the rates. The DSO theory has been proven to be 
quantitatively valid down to temperatures 4kgT ~ y 0 
in the single electron transistor [23]. Upon decreasing 
of the temperature below 70 /4, a quantitative descrip¬ 
tion of the transition rate r“ c would require to calculate 
E beyond the lowest order in j 0 . In the regime where 
temperature and coupling are of comparable magnitude, 
the width and position of the Coulomb blockade peaks 
in a gate trace are strongly influenced by the TDOS and, 
more precisely, by the self energy E. The role of Re(E) is 
to influence the positions of the Coulomb blockade peaks: 
In the rate for the transition a to b , the real part appears 
next to the energy difference E b of the transition in the 
denominator. Hence, due to this contribution the res¬ 
onant level is shifted depending on the configuration of 
the leads. 


D. Renormalization of excited states 
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FIG. 7. Example of possible charge fluctuations for a final 
state (|AT), left, shaded gray) with one extra electron and 
an initial state with zero electrons in the shell n (| 0 ), right, 
shaded gray). This set corresponds to one specific diagram 
of the type shown in Fig. IT). States that can be reached 
by in-tunneling of an electron are shown on top, states that 
can be reached by out-tunneling of an electron are shown on 
the bottom. Dashed frames highlight resonant (E°j b = 0) 
charge fluctuations. Above and below the level schemes, the 
energy difference between the virtual state and the state on 
the other contour is given: the energies of the states accessible 
from the initial (final) state are compared to the energy of 
the final (initial) state on resonance ( E^ = eaV g ). Note 
that the electron number of the states that can be reached 
by in-tunneling on the left and the number of electrons in the 
initial state on the right differ by two. The same situation 
occurs for the final state and the out-tunneling states on the 
right. The energy differences for this class of fluctuations is 
of the order of E c . A comparison of the electron number of 
the final state with the in-tunneling states on the left and the 
initial state with the out-tunneling states on the right yields 
a difference of zero. These fluctuations have comparably low 
energy cost. 


processes. Although the framework in Ref. [32| is differ¬ 
ent, the evaluation of Re(E ba ) is similar to that in our 
model. The condition for a resonance for a transition 
between states a and b is given by 

e ± eV b /2 + eaV g - E b + Re(E ba ) = 0, (14) 


In the stability diagram in Fig. [2] we observe an asym¬ 
metry in the spacing of lines associated with excited 
states connected to one charging state, as drawn schemat¬ 
ically in Fig. [s| The line 0 —> V meets the diamond at 
bias voltage I4i- Measured along the bias voltage axis, 
this value is larger than the energy difference Vb 2 asso¬ 
ciated with the line 2 —>• V on the right. A similar be¬ 
havior has been discussed previously for the co-tunneling 
regime [32]. As noted by these authors, the asymmetry 
can not be explained within the sequential tunneling pic¬ 
ture but can be attributed to the renormalization of the 
excitation energies E b in Eq. ([ 9 ]) due to virtual tunneling 


where e is the energy of the tunneling electron with re¬ 
spect to the chemical potential of the unbiased contact 
fi 0 . Note that this condition can be fulfilled for differ¬ 
ent transitions at the same time, a situation that occurs 
at any point where two lines in a stability diagram in¬ 
tersect. In order to interpret the observed shift of the 
excited state line in the differential conductance data in 
Fig.0 it is illuminating to study the contribution from 
Re(E) at points (k^i, I4i) and (V g 2 , Vb 2 ) marked by a dot 
and a circle, respectively, in Fig. [8] We consider an exem¬ 
plary set of states 0 = |0;n), 1 = |A'f h), li = |[Af];n), 
I 2 = |A'f) AT4-, (AT); n) and 2 = |Af, Ki;n). A similar 




































FIG. 8. Schematic drawing of the conductance lines in the 
vicinity of the charging state with 4n + 1 electrons in Fig. [2] 
The first visible excitation is shifted upwards on the left and 
downwards on the right side of one charging diamond by 
e(Vbi — Vb2)/2 = —Si. The corresponding energies in Fig. [2] 
are eVbi/2 ~ 2meV and e 14 , 2/2 — 1.4 meV. For our analysis 
we choose bias and gate voltages close to the filled dot for the 
first transition 0 —> 1 ' and to the empty circle for the second 
transition 2 —» lb 


analysis can be carried out for other states with 4n + 1 
and 4n + 2 electrons. The quantum numbers in round 
brackets denote a missing electron of shell n—l whereas 
the square brackets indicate a state of shell n + 1. For 
each of the highlighted points in Fig. [8j two conditions 
in the form of Eq. (14) can be given. Subtracting them 
pairwise we are left with 


eV hl - E \' + [Re(E 1,, °) - Re(E 1,0 )] = 0, (15) 
eVb 2 - El' + [Re(E 2,1 ) - Re(E 2 ’ 1 ')] = 0, (16) 


where the self energy contributions depend on bias and 
gate voltage. To lowest order in 70 we analyze the dif¬ 
ferences in Re(E) using eVbi/2 = E}' and okW ^ 2 = 
Eyjf ± eVb /2 at e = 0. In order to calculate Re(E) 
we have to analyze the contributions from all acces¬ 
sible states in Eq. (10). In principle there are arbi¬ 
trarily many states that can be reached by a charge 
fluctuation. However, we assert that the available en¬ 
ergy interval for charge fluctuation processes is given by 
max(eVb,F 0 ,3 — 4keT) and contributions beyond this 
scale are suppressed. Numerical results using a larger 
bandwidth can be found in Sec. [A] of the appendix. 


For our considerations we assume that the spin orbit 
coupling of our CNT quantum dot is small, i.e., Ago < 
max(keT,r). Otherwise we would expect to see a two¬ 
fold symmetry in the spacing of the excited state lines 
in the stability diagram in Fig. [2] The other important 
scales - charging energy, shell spacing and linewidth - are 
related in the way E c > e 0 max(kBT,70). Within 
this choice of parameters the difference of the self energy 
corrections for the resonant transition can be calculated 


by (15) —(16)= 0, i.e., 


61 = 


Re(E 10 ) - Re(E 10 ) 


Re(E 2i ) -Re(E 21 ') 


— 7o 1 + 2 k s —Kd + Kf — kj.} ’I , b(£o/2) 


( 17 ) 


where we used the abbreviation = Re[’F°(l/2 + 

ie/27rkBT)] and a bar denotes a summation over indices, 
e.g., Ri = Kicr- A detailed derivation of these quanti¬ 
ties is given in the appendix, Sec.[E] Similar calculations 
are performed for the excited states in the n + 2 and n + 3 
diamonds, yielding 

62 ^ 7 o {its-Kd + K s .|.-Kdf} 'Fr( £ o/ 2 ), 

S 3 ~ 70 {1 + K s -2S, d + K4.-%}T , B(e 0 /2), 


where the states with three electrons are chosen to be 
electron-hole symmetric with respect to the state with 
one electron. Note that for the case of symmetric cou¬ 
plings the shifts reflect the electron-hole symmetry of the 
system while a choice of a 7^ 0 (Eq. <§) breaks this sym¬ 
metry. For highly asymmetric couplings |a| ~ 1 the shifts 
are comparable to those in Ref. [32]. Note that the ef¬ 
fective change of the resonance with respect to the en¬ 
ergy difference has a negative sign (compare Eq. ©)• 
The resonance marked by the left arrow in Fig. [2] is sit¬ 
uated above the resonance marked by the right arrow. 
The experimental data thus corresponds to a negative 
shift. We therefore assume an asymmetric coupling to 
the leads with a dominant coupling to the drain contact, 
i.e., k s < K d , — 1 < a < 0. Using the parameters from a 
fit to the data in Sec.[V] i.e., a = —0.7 and eo = 1.4meV 
we obtain <5i ps —0.2meV and 62 ~ —O.lrneV. Compared 
to the shifts in the experimental data, these values are 
too small by a factor of 2-3. We expect that additional 
states may contribute to the charge fluctuations that are 
not considered within this approximation. 


E. Tunneling magneto-resistance 


Corrections to the conductance peak width are given 
by Im(E). Because Re(E) and Im(E) both depend on 
the different magnetic properties of the source and drain 
leads as well as on the dot’s configuration, the resulting 
impact on the TMR is quite intricate. Thus we analyze 
the contributions to the self energy in the light of different 
configuration of the lead’s polarizations. We focus on the 
last resonance, i.e., the transitions |0, n+1) {|(crr), n+ 

1)} where the TMR graph in Fig. [4] exhibits a double dip 
like structure. The back-gate voltage is tuned such that 

e + eaV e -E° {Ta) +Re( E 0 ’^) = 0, 


and the quantum numbers in round brackets (rer) denote 
a missing electron of shell n + 1. At lowest order in 
the tunnel coupling 70 we approximate eaV g = E® Tcr ^ 

when we calculate Re(E°d T<T )). From Eq. (12) we list the 
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imaginary part of the self energy for this transition, i.e., 

Im(S°’^))=7r 7n ^ 

l 

| K l<T(c)ff{E( T <r) ~ e ) + y K lcr(c')fi (e — ^(tct)) 

v cec+ c'ec~ 

+ K l<r(c)ff( e ~ -S'c) + y~l K la(c')fl (-E’c'^ e )|- 

cec+ , c'ecr . J 

(TO-) 


The magnitude of the energy difference of the virtual 
state with respect to the state on the other contour deter¬ 
mines whether a possible charge fluctuation contributes 
to the renormalization of the self energy or not: a contri¬ 
bution f+(E c — e), e.g., is exponentially suppressed in 
the vicinity of the resonance. 

Therefore, knowing the arguments in the step func¬ 
tions / , we can simplify the result significantly. Close 
to the resonance where |e| < max(keT, 70 ), the fluctua¬ 
tions with an energy cost of the charging energy E c or of 
the shell spacing e 0 , e.g., the states that can be reached 
by out-tunneling from the state |(tct)) can be neglected. 
Focusing on the resonant contributions, we are left with 

—- - H ( e ) + Y, ( e - E t)^> } • 

0 1 { J 

(18) 



It is clear from this result that the broadening of the 
TDOS peak does depend on the lead configuration {ni a }. 
Let the majority spins be polarized such that a = +1 
in the layout with parallel lead polarization. The sum 
over the leads is then given by J2i K fa = (1 + oP)j 4 
and Yli = (1 + nPa)/4 for parallel and anti-parallel 
polarizations, respectively. Let us first consider the case 
of zero effective Zeeman splitting, i.e., E ? = = 0. 

The difference of Im(E) for the two configurations then 
reads 


Irn 


-^0,(tct) 


- E 


ap 


= S lm = 7T^aP(l-a)f + (e). 


(19) 


Note that the validity of this result depends on the ra¬ 
tio of linewidth and level spacing, namely that 7 o -C eo 
such that only the selected small set of charge fluctua¬ 
tions contribute. The sign of the difference in Eq. (191 is 
determined by cr, a result which is intuitively clear since 
the sum over the couplings will be greater for the spin-up 
transition (a = 1 ) in the parallel case and for the spin- 
down transition in the anti-parallel one (a = — 1 ), as 
shown schematically in Fig. [9](a). For zero energy split¬ 
ting E^j we would expect a broadening of the peak as¬ 
sociated with the transition 0 (f) for the parallel con¬ 

figuration and a broadening of the peak in G ap for the 
transition 0 (j.). Note, however, that the second effect 

will not be visible since the TMR ratio will be dominated 



FIG. 9. (Color online) The influence of Im(E) on the TMR. 

a) Large gray arrows symbolize the majority spin in the left 
or right contact. The contributions to the self energy for one 
spin species are summed for each configuration of polarized 
leads (parallel on the left, anti-parallel on the right) as indi¬ 
cated by the dashed frames. Weak (strong) coupling to the 
dot (blue ellipse) is given by thin (thick) arrows. Note that for 
the spin down species the sum over the leads yields a greater 
contribution in the configuration with anti-parallel polariza¬ 
tion (as indicated by the signs between the dashed frames). 

b) On the left, we depict schematically the conductance peaks 
for one resonance in both parallel and anti-parallel configura¬ 
tions and the resulting TMR (right). The broadening of G p 
is typically larger than for G ap in the absence of stray fields. 

c) Due to a magnetic stray field, the contribution to Im(E) in 
the parallel case can be reduced, giving rise to a double dip 
structure in the TMR. 


by the spin up transition. Hence, we will observe a TMR 
signal as depicted in Fig. §>)• 

Now let us assume a non-zero effective Zeeman split¬ 
ting E] = £7 — Ei = gfrBft.p/ a p of states with quan¬ 
tum numbers a =f / 4- This splitting also depends on 
the magnetization state p (parallel) or ap (anti-parallel) 
of the contact electrodes. The energy difference is ex¬ 
pressed in terms of the effective magnetic fields g^h p 
and g/ie^ap- We assume that this field is non-zero for 
both polarizations. Im(E) as well as the TMR are very 
sensitive to the choice of the shifts, the couplings and 
the polarization. The mechanism we want to discuss can 
be observed for different parameter regimes, but for the 
sake of the argument it is sufficient to present one pos¬ 
sible set that we deduce from the experiment and the 
line of reasoning that goes with it. In the last part of 
Sec. |III D| we argue that couplings k s < or, similarly, 
0 > a > — 1 are needed to explain the shift of the ex¬ 
cited state lines in Fig. [2j Furthermore we point out 
that the peaks in conductance in Fig. [I] are descend¬ 
ing in height as we fill the shell. In our model the 
drain lead switches polarization upon interaction with 
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external magnetic field while the density of states in the 
weakly coupled source contact remains unaltered. Given 
that the spin transport is more sensitive to the bottle¬ 
neck (source) contact, it is plausible to assume that the 
shifts are such that the majority spins tunnel first on 
the quantum dot, namely spin up electrons in both con¬ 
figurations. These considerations favor a choice of neg¬ 
ative shifts g/r.Bh a p, gp^h-p < — kpT. The second pair 
of resonances is then dominated by spin down electrons 
and the respective contributions f~(e + ggB^-p/ap) in 
Eq. ( fl8| ) are suppressed. Conversely, for spin up electrons 
/“(e — ggB^p/ap) = 1- In the resonant case, |e| < keT, 
the imaginary part of the self-energy for the |0 ) ^ |(ct)) 
then reads 

5 Im ^-7r|(l-a)P(l-a/-(e)). (20) 


The magnitude of the relative broadening of the peak re¬ 
lated to the transition of a spin down electron in G ap is 
thus increased for higher polarization and a —> — 1. Al¬ 
though this estimate is only valid in the direct vicinity 
of the resonance, it describes the situation qualitatively 
as can be seen in Fig. [lO] We show conductance and 
TMR nearby the resonance |0,?z + 1) ^ {|(crr),n + 1)} 
for fields ggs/ip = —40 geV and ggs/iap = — 80 geV. In 
the panels on the left side, the polarization is varied keep¬ 
ing a = —0.8 fixed. We see that the right shoulder in the 
TMR curve (c) is lifted upwards with increasing polar¬ 
ization. On the right panels in Fig. [To] we increase the 
coupling to the source contact which is proportional to 
a. While the conductance is decreased for asymmetric 
choices of a in both configurations (see (d) and (e)), the 
magnitude of the peak in G ap is not symmetric with re¬ 
spect to the coupling to source and drain. The TMR in 
Fig. [Hjtc) can be related to Eq. (20): the shoulders for 


a = 0.8 turn into dips approaching a = —0.8. Please 
keep in mind that this discussion is simplified since we 
do not account for the fact that the relative position of 
the peaks changes, too, as we vary the parameters a and 
P (compare Re(S) and Im(E) plotted in Fig. 12 in the 
appendix Sec. [A|. 


IV. COMPARISON 


a. Con ductance in the experiment and in the model 
In Fig.lllla) (blue circles) we show the conductance Gp ast 
obtained at B = 0 performing a fast measurement, i.e., 
sweeping the gate voltage V g at zero bias voltage, see 
Sec. |IIB| Note that it provides only conductance data 
for the parallel configuration (compare Fig. [3]) . The data 
from this measurement yields conductance peaks that fit 
to Lorentzian curves with an average FWHM of 0.3 meV. 
Adapting our model parameters to the data of Gp ast , 
we obtain the continuous lines in Fig. ll'a,b,d). The 
cond ucta nce data from the slow measurement (compare 


Sec. IIB) for the two configurations, Gp low and G® 1 "", 
are shown in Fig. |TT|(a,b) (green crosses). The shape of 






FIG. 10. (Color online) Conductance and TMR calculations 
in the vicinity of the resonance |0,n + 1} {|(crr),n 4- 1)} 

for different polarizations P (panels (a)-(c), a = —0.8) and 
coupling asymmetry a (panels (d)-(f), P = 0.4) applied 
in the parallel configuration for effective Zeeman splitting 
gUshp = —40geV and ggu/iap = —80 geV. (a),(b): increas¬ 
ing the polarization reduces the peak width and height of both 
G p and G ap . (c): In the TMR curve, the shoulder on the left 
at P = 0.2 is shifted to the right for P = 0.6. (d),(e): The 
coupling asymmetry a^0 diminishes the peak heights of the 
conductance for both configurations of the leads. Note that 
in the anti-parallel case shown in (e) the symmetry between 
the contacts is broken and the peak height is sensitive to the 
variation of the dominating coupling, (f): The TMR curve 
exhibits a double dip feature for values — 1 < a < 0. It is 
transformed to a double peak for 0 < a < 1. All plots are 
calculated at a temperature corresponding to 40 geV and a 
coupling 7 q = 160 geV. 


the conductance peaks turns out to be non-Lorentzian, 
with the peak height in the conductance data limited to 
~ 0.1e 2 /h. While the flanks of the peaks match for the 
first three resonances in the data from the slow and from 
the fast measurement EDI. the maximum conductance 
values deviate by a factor of three. So far no full expla¬ 
nation for the suppression of the peak conductance was 
found. 

b. Model parameters A bare coupling of 70 = 

80geV is found to optimize the fit to Gp ast . The ther¬ 
mal energy is chosen as JcbT = 40 geV (460 mK), close to 
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FIG. 11. (Color online) Conductance at zero bias as a function of gate voltage V g plotted for (a) parallel and (b) anti-parallel 
polarization of the leads. In (a), a gate trace (Gp ast (Vjr), blue circles) is shown together with conductance obtained during 
TMR measurements Gp low (B,V g ) (green crosses, see also Fig.^J), and the calculated conductance for parallel lead polarization 
(continuous line, black) at kaT = 40 ^teV, eo = 1.4 meV, E c = 6.1 meV, a = —0.7 , P = 0.4, g/m/iap = —0.16 meV and 
gfishp = — 0.12meV. In the vicinity of the rightmost resonance, Gp ast shows a high noise level (compare also Fig. [ 2 J). (b) The 
conductance data measured for anti-parallel polarization of the contacts Gap W (B, Vg) (green crosses) is compared to the model 
output (continuous line, black) for the same parameters as in (a), (c) Experimental TMR data calculated from Gp low (a) and 
Gap W (b) (also shown in Fig. (d) TMR obtained from the model conductance (continuous lines in (a) and (b)). 


the base temperature (300 mK). For the quantum dot 
parameters we set E c = 6.1 meV and a shell spacing 
e 0 = 1.4 meV as inferred from Sec. |TT] The shell num¬ 
ber n ~ 40 is estimated from the distance to the band- 
gap. We assume asymmetric contacts with a = —0.7 
and polarization P = 0.4. For the calculation of the 
charge fluctuations we include all states within an en¬ 
ergy interval of 3eo (see Sec. [A] in the appendix). The 
effective Zeeman shifts for the model output in Fig. [Tl] 
are gg b/i p = — 0.12meV and g/-iBh ap = — 0.16 meV. 


c. Discussion If only features of the leads density of 
states at the Fermi energy are included, compare Eq. 
the DSO preserves particle-hole symmetry by construc¬ 
tion 22- To break this symmetry, a Stoner-shift of the 
majority band with respect to the minority band should 
be included [25] , whose effect is analogous to that of an 
effective Zeeman field [21]. Such effective fields have also 
been used to model the effects of coherent reflections 
at the magnetic interfaces in double barrier systems [7|. 
Since the data in Fig. 11 'a-c) does not reflect particle hole 
symmetry, we use effective Zeeman splittings to break the 
particle-hole symmetry and reproduce the observed mag¬ 
nitude of the TMR effect. The splittings are of similar 
magnitude as those used in Ref. HQ] (g/i b/i p = 0.25 meV 
and ggsh ap = 0.05 meV) to explain the experimental 


TMR data of Ref. 0 

In case of non-zero spin-orbit coupling sd sa , we 
would expect a splitting of the excited state lines in the 
stability diagram in Fig. [2] This is not resolved in our 
experimental data. For simplicity we therefore here as¬ 
sume Ago = 0. Model calculations with non-zero spin 
orbit coupling can be found in the appendix, Sec. [C] 

From the conductance traces calculated within our 
model, Fig. |TT)(a,b) (continuous lines), the TMR, 
Fig. [Tljd) , is obtained. The data and the model calcu¬ 
lation agree in the decay of the TMR amplitude within 
a sequence of four charging states including the “double 
dip” feature in the last two resonances at V g = 8.17 V and 
V g = 8.19 V. This indicates that the sequence in Fig. [IT] 
represents one shell, i.e., charging states 4n+l to 4(n+ly 
We note that in the model output the last resonance is 
dominated by a peak while the dips are more prominent 
in the experimental data. 

In the vicinity of all conductance peaks (at V g = 
8.13 V, V g = 8.15 V, V g = 8.17 V and Vg = 8.19 V) an 
additional small shoulder around TMR = 0 occurs in the 
data of Fig. E3 c). These shoulders are likely related to 
the aforementioned suppression of the peak conductance 
in the slow measurement (see Fig. ©a ,b)). We recall 
that the TMR is calculated from the ratio G p /G ap (com- 
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pare also Fig. [5] and Fig. [9]): in the regions where the 
peaks are cut off, the ratio Gp low /G a ° w is smaller than it 
is in the same region in the model output, where steep 
peak flanks lead to a larger ratio G p /G ap . 


V. SUMMARY 

The tunneling magneto-resistance of a carbon- 
nanotube based quantum dot with ferromagnetic leads 
has been explored both experimentally and theoretically. 
The experimental data shows a distinct variation of the 
tunneling magneto-resistance (TMR) lineshapes within a 
single quadruplet of charging states. 

To model the data we apply the dressed second-order 
(DSO) framework based on the reduced density matrix 
formalism. This theory accounts for charge fluctuations 
between the quantum dot and the ferromagnetic con¬ 
tacts. Thereby, it goes beyond the sequential tunnel¬ 
ing approximation which can only account for a positive 
and gate-independent TMR. When the charge fluctua¬ 
tion processes are summed to all orders in the coupling 
to the leads according to the DSO scheme, they yield 
tunneling rates where the Lamb shift and the broaden¬ 
ing of the resonances are given by the real and imaginary 
parts of the self energy, respectively. This is a nontrivial 
result which yields the tunneling rates for an interact¬ 
ing quantum dot in the intermediate parameter regime 
E c ke;T ~ r depending on the polarization of the 
contacts. 

We explicitly compare the DSO self energy for different 
contact magnetizations and show that the DSO model¬ 
ing can account both for the renormalization of excited 
states and the specific structures observed in the TMR 
gate dependence. A comparison of the TMR obtained 
from the model and from the experimental data shows a 
qualitative agreement. 


ACKNOWLEDGMENTS 

We gratefully acknowledge discussions with Johannes 
Kern, Davide Mantelli and Daniel Schmid. This work 
was funded by the Deutsche Forschungsgemeinschaft 
(DFG) via GRK 1570, SFB 689 and Emmy Noether 
project Hu 1808-1. 


Wfluc — 70 . e o . 2e 0 3eo 

Re(S p )/ 70 ~ 6k * T : Re(£ ap )/ 7 o 


Re(S p )/ 70 ~ 6k * T : Re(£ ap )/7o 



- 2-1012 e / f0 - 2-101 2 


FIG. 12. (Color online) Re(E) (a,b) and Im(S) (c,d) for both 
lead configurations as a function of energy t in units of the 
shell spacing eo- Different lines are plotted for bandwidth 
Wfl U c = 70 (green, dotted) to 3eo (red, continuous) in steps of 
eo• In the vicinity of a few keT around the resonance (e = 0, 
gray region) the difference between the graphs for the real part 
(a,b) is small and for the imaginary part (c,d), it is vanishing. 


available energy scale in the system. Nevertheless it is in¬ 
teresting to see how the outcome is affected by increasing 
the bandwidth and allowing excited states of the neigh¬ 
boring shells to contribute to the charge fluctuation chan¬ 
nels. In terms of an effective energy shift in a multi-level 
quantum dot the renormalization due to excited states 
was also discussed in Ref. EH To illustrate the effect 
of such a modification we plot the real and imaginary 
parts of the self energy T, in the vicinity of the transition 
|(K i),n) |-,n + 1) for different sets of charge fluc¬ 

tuations within energy ranges of 70, ecu 2eo and 3eo in 
Fig. [l2j We clearly see that the fluctuations from higher 
shells manifest themselves in additional features in the 
curves for Re(£), Fig. 12 (a,b), and Im(E), (c,d). Note, 
however, that the zero-bias conductance in our system 
is only sensitive to a small vicinity of a few kgT around 
the resonance. Within this range the high energy con¬ 
tributions do not change the picture substantially. The 
analysis of the imaginary part in Sec. |III E| is thus exact 
at the level of the self energy since the Fermi functions 
in the imaginary part suppress contributions from other 
shells. 


Appendix A: Contribution of other excited states to 
the renormalization of the self energy 


Appendix B: Calculation of Re(E) 


Whenwe discuss theeffect of the charge fluctuations In this section we perform the calculation of Re(E 1 ' 0 ) - 
Sec. HID and Sec. HIE of the mam text, we always twvi<L x x Ijno 


in Sec. HID and Sec. HIE of the main text, we always 
focus on the most resonant transitions (see Fig. [?| that 
are energetically favorable, i.e., on transitions in Eq. fl2| 

c' /b ' -* 

with an energy difference E a j c , of the order of the effec¬ 
tive line-width or below. At zero bias this is the largest 


Re(£ 10 ) as part of the quantity Ji introduced in Sec. 
of the main text. To this extent we analyze the renormal¬ 
ization of the energy difference E\ due to charge fluctu¬ 
ations to and from states 0 = |0;n), 1 = |/\f;n) and 
1' = | [Af]; n) in more detail. We recall that the real part 
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of the self energy related to a charge fluctuation to state 
d has the form (see. Eq. e§) 

-E-rf ( * / “ ) W®&(w±p(<2-0). 


where we have to replace 6 =l , ,a = 0 or 6 =l and a = 0 , 
respectively. Note that the contribution oc v I ;(0 \l-E) in 
Eq. (12) does not appear explicitly since it cancels in the 


difference of the shifts. Next, we have to find all states d 
that contribute within our resonant approximation. We 
can immediately discard states that can be reached by in¬ 
tunneling from b and by out-tunneling from a , since their 
energy differences E c a j ^ are of the order of the charg¬ 
ing energy and thus beyond our charge fluctuation band¬ 
width of W e = max(eVb,kBT, 7 o) = eo/2. We are left 
with states that can be reached by in-tunneling into state 
a and by out-tunneling from state b. Let us discuss one 
example for the state 1'. There is one electron in the shell 
n +1 (denoted by the brackets [...] in the state ket) which 
can tunnel out and we are left with a state |-, n). Actually 
this state is identical to the state 0 on the other contour, 
thus Eq =0 = 0. We can now evaluate the argument of 
the digamma function, i.e., /q — Eq + e, for e = 0. Since 
Hs/q = ±eo/2 and thus |/q| < W e , we have to sum over 
both leads. The total contribution from fluctuations to 
d = 0 is thus —70 K it 4-^(e 0 /2). The other states that 

can be reached by out-tunneling, e.g., |(A"f), [Kf], n), 
yield energy differences of at least 3/2e 0 > W e . Using 
similar arguments we can collect all relevant contribu¬ 
tions to the difference Re(E 1 °) — Re(E 10 ). In a graphical 
representation, this can be visualized as 


in to 0 


Re(E 10 ) - Re(E 10 ) = 


-J330UT 

out from 1 ' 


JjD000" 


out from 1 



0 

QO© 


bed 

U £( X~ 7 



= 2k s 4-0 (e 0 /2) 

where one set of four boxes symbolizes one shell and 
we use E c ! b . as a label. Fluctuations that cancel are 

a/c' 

crossed out. Note that for excited states with an energy 
difference E c a j b = ±eo we add only the contribution 

from the source (drain) contact where |/q — E° a j b ,\ < W e . 
Similarly we find 

Re(S 21 ) - Re(E 21 ') = (1 + R d - % + ^)4-^(e 0 /2), 


2 

TMR 

1 


0 


8.13 8.15 8.17 8.19 

gate voltage (V) 

FIG. 13. TMR as a function of gate voltage for orbital po¬ 
larization P OT b = 0.6, orbital shifts £r 0 rb/Torbhap b = —80 peV 
and <?orbMorbhp = — 40 peV, and Aso = 0.1 meV at k B T = 
40 peV. The other parameters are identical to the ones used 
in Fig. [TT] 



which leaves us with di from Eq. (|17[). 


Appendix C: Spin-orbit coupling and valley 
polarization 


In Sec. |III A| we discussed the possibility to include 
spin-orbit interaction effects, as they have been reported 
to play a prominent role in carbon nanotubes mm- 
However, we did not add it in the comparison to the ex¬ 
perimental data since they could not be resolved in the 
transport spectrum (Fig. [ 2 ]). Nevertheless, values of the 
order of Ago ~ 100 fj,eV would still be consistent with 
the experimental data. Introducing a finite Ago a priori 
does not affect the TMR as the Kramers pairs are spin de¬ 
generate pairs with anti-parallel and parallel alignment 
of spin and valley magnetic moments. Yet it has been 
argued that the two valleys of a CNT can couple differ¬ 
ently to the leads [53] , If the valley quantum number 
is conserved upon tunneling, the mechanism can be un¬ 
derstood in terms of a valley polarization. A possible 
tunneling Hamiltonian that describes this situation can 
be written as 

Ht = E dj tcrr Cgka- “b h. C., (Cl) 

Zknerr 


with a valley dependent coupling T;k, lcrr and an opera¬ 
tor cjkrtr that describes the electrons in the leads (that 
are also part of the CNT). Including a valley polariza¬ 
tion in turn also renders the TMR sensitive to mag¬ 
netic stray fields g 0 rbMorbfrp rb and 5 orbAtorb/iap b alon g the 
tube axis. The orbital magnetic moments gorbMorb are 
considered to be larger then /in by one order of mag¬ 
nitude 23- I n Fig- 13 we present a TMR calculation 
for Ago = 100/reV, orbital polarization P orb = 0.6 and 
stray fields 5 orbMorb^ap b = -80^eV and g 0 rbf«orb^p rb = 
—40 /reV again combined with the experimental data. 
The spin-dependent shifts are assumed to be negligible 
in this setup. We see that the agreement with the exper¬ 
imental data improved slightly in Fig. 13 at the expense 
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of additional free parameters. It is, however, outside the scope of this paper to discuss the effect of spin-orbit cou¬ 
pling and the valley polarization in more detail. 
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